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Abstract 

The classical Yang-Mills equation of motion is numerically investigated in 
the Lorentz gauge for a SU(2) gauge group. The color-electric field of two 
point-like charges is studied in the "empty" vacuum and in a state with an 
instanton present. The major effect for a fixed orientation of the instanton is 
that the color-electric field lines are expelled or attracted from the instanton 
region depending on the orientation of the instanton. If over the orientations 
of the instanton is averaged, this effect drops out. In this case of a random 
instanton orientation, we find that the external color-electric field is expelled 
from the instanton core. The origin of this effect is discussed. 
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1 Introduction 



One of the most challenging problems in hadron physics nowadays is to understand 
the confinement of quarks in quantum-chromodynamics (QCD), which is the theory 
of strong interactions. From large scale numerical investigations of lattice Yang- 
Mills theory |l], Q], we know that QCD yields confinement, since Wilson's area law is 
satisfied. This implies that a linear rising confinement potential is formed between 
static color sources. Unfortunately, the numerical efforts have not yet revealed 
the basic mechanism of confinement. The knowledge of the qualitative mechanism 
would help to construct effective hadron theories the predictive power of which is 
not limited by unphysical quark thresholds ||. Only effective quark models which 
heuristically incorporate the confinement are so far available. I would like to mention 
the Global Color Models @ which realize confinement either by avoiding the pole at 
the real axis of the quark propagator or by infra-red slavery. In || it was proposed 
that quark confinement is due to random quark interactions induced by random 
gluonic background fields. 

A promising idea to understand the feature of quark confinement from first prin- 
ciples was given by 't Hooft and Mandelstam. They propose that the so-called 
Abelian gauges are suited to understand the basic mechanism of confinement which 
is obscured in other gauges [El. In these Abelian gauges magnetic monopoles emerge 
because of a residual U(l) gauge degree of freedom which is unfixed 0. If these 
monopoles condense due to the Yang-Mills dynamics, a dual Meissner effect takes 
place implying that the color-electric flux is expelled from the vacuum. This scenario 
would naturally explain the linear rising potential between sources of color-electric 
flux. Subsequently, heuristic models were developed which explore the formation of 
an color-electric flux tube [§]. However, the dynamics which leads to the condensa- 
tion of the monopoles is not understood so far. 

The original idea of 't Hooft and Mandelstam has gained further support by a 
recent success || of Seiberg and Witten. They showed that the monopoles of N = 2 
(non-confining) super-symmetric Yang-Mills theory start to condense, if the theory 
is explicitly broken down to N = 1 (confining) SUSY. One should, however, keep in 
mind that this scenario in super-symmetric Yang-Mills theory is quite different from 
that in standard QCD. Whereas in the previous case magnetic monopoles are present 
in the particle spectrum, magnetic monopoles are induced by gauge transformations 
in the latter case. 

The quasi-particles of standard Yang-Mills theory which show up as a solution of 



the classical Yang-Mills equation of motion are instantons [10|. An implicit scheme 



to construct all instanton solutions was provided in [11]. Novel explicit instanton 



configurations for a SU(N), N > 3 can be found in ]P2]|. Instantons are generally 



believed to play an important role in the QCD ground state. They may possibly 
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trigger the spontaneous breaking of chiral symmetry |L3] and offer an explanation 
of the £7a(1) problem |TJ . 

Numerical investigations, based on a classical instanton interaction, show that 
instantons are strongly correlated in the Yang-Mills ground state implying that 
they occur as a liquid [O rather than in a dilute gas phase 1TH] . Including quantum 



fluctuations, it was discovered that instantons possess a medium range attractive 
interaction which is solely due to the instability of the perturbative vacuum flT? . 



The strong-coupling expansion within the field strength approach to Yang-Mills 
theory [|18[ suggests that the interaction due to quantum effects is strong enough to 
induce a condensation of instantons [ 19(| . Whether the instantons exist in a crystal 
type structure or in a strongly correlated liquid is still beyond the scope of the 
present approaches. 

The question whether instantons are important for the dual Meissner effect has 

In order to address this question 



gained recent attention [2C, 21 



H li- 



the 

profile of an instanton in Abelian gauges were investigated in continuum Yang-Mills 
theory |2(| |2T|] as well as in the lattice version [22, |23 



In 20 



a monopole world line 
was observed which penetrates the center of the instanton. In contrast, the lattice 
have reported closed loops of monopoles which encircle the center 



simulations 22 



of the instanton. The length of the monopole loop was investigated throughout the 
deconfinement transition. The discrepancy between the lattice and the continuum 
approach was clarified by Brower at al. pl| |. They showed that the configurations 
which minimize the gauge fixing functional indeed lead to closed monopole world 
lines. The conclusion at hand is that large monopole loops play an important role 
for confinement. Recently, lattice calculations reported an enhanced probability for 
monopoles inside the instanton |24| indicating that instantons play a role in the 
context of the dual Meissner effect. 

In this paper, we will directly focus on this role of the instantons. For this 
purpose, we will numerically solve the SU(2) Yang-Mills equation of motion with 
two static color sources present. The color-electric flux stemming from the color 
sources will be investigated in the "empty" vacuum and in a vacuum where an 
instanton is present. Whether the external color-electric field lines are expelled or 
attracted depends on the orientation of the instanton relative to the charges. If over 
the orientation of the instanton is averaged, the picture changes qualitatively. We 
will provide evidence that in the latter case the external field is repulsed from the 
instanton core. 

The paper is organized as follows: in the following section, we briefly review 
some aspects concerning the relation of instantons and confinement. In section |3|, 
the Meissner effect in solid state physics is addressed for a later comparison with 
the situation in SU(2) Yang-Mills theory. The numerical approach to solve the clas- 
sical Yang-Mills equation of motion is introduced in section It is shown that 
this approach nicely reproduces an instanton configuration. In subsection 4.3, two 
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point-like color-electric charges are added. The electric field lines stemming from 
these charges are investigated in an "empty" vacuum (subsection fO| ) and in a state 
with an instanton present (subsection fOj). The case of a random instanton orien- 
tation is presented in section Conclusions as well as an argument to understand 
qualitatively the numerical results of the latter section are left to the final section. 



2 Instantons imply confinement? 

In this section, we briefly review some results in the literature which address the 
question how relevant instantons are for the confinement mechanism. 

This question has a long history |25|, ^6|, 27]]. A semi-classical evaluation of the 



Yang-Mills partition function by Aharonov et al. p5] has indicated that instantons 
support the existence of infinite long flux tubes (fluxon). These strings with non- 
vanishing flux appear as fluctuations on top of the instanton background. It was 
shown that an ensemble of these random- walking "fluxons" gives rise to Wilson's 
area law. In this scenario, the impact of the instantons on confinement is secondary. 
Considering fluctuations around an instanton, 't Hooft showed that large size 



instantons receive a large weight in the Yang-Mills partition function [14]]. Expec- 
tation values of condensates are even ill-defined due to the contribution of these 
large size instantons. Subsequently, it was argued by Shuryak that the presence of 
quarks and gluons cure the large size instanton problem |[26|| . He showed that the 
presence of a non- vanishing quark and/or gluon density provides a natural cutoff 
on the instanton radius. As a consequence, the presence of quarks and gluons yield 
an increase in the vacuum energy density implying that they should be expelled 
from an instanton dominated vacuum. This picture naturally supports the MIT bag 
model of hadronic matter. Note that we will focus at the classical level when we 
will elaborate below the response of the instantons on the external charges. 

In their pioneering work, Polikarpov and Veselov provided evidence that the 
instantons play a minor role for confinement f2~7fl . They studied the emergence of in- 
stantons in SU(2) lattice Yang-Mills theory, when a equilibrium lattice configuration 
is cooled down. They extracted the ratio \ of the string tension of the configuration 
after cooling with the one of the lattice equilibrium configuration as function of the 
topological charge present on the lattice. One crucial observation is that \ increases 
with increasing topological charge Q. For instance, x reaches 20% for a sufficiently 
large Wilson loop and for Q = 3. However, large topological charges rarely appear 
on the lattice. The value of x averaged over the instanton vacua which appear by 
freezing out equilibrium configurations is of order 5%. 

The outcome of the latter considerations is twofold: an instanton medium which 
emerges from the freezing of a lattice equilibrium configuration cannot give rise to 
the full string tension. Fluctuations on top of the instanton medium seem to play 
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an important role for confinement. Recent lattice simulations indicate that these 
fluctuations might be magnetic strings, which are attached to the instanton [EM . 



Secondly, lattice equilibrium configurations are dominated by fluctuations rather 
than by the classical solutions of the Yang-Mills equation of motion. The later fact 



is confirmed by recent studies in p3]| . The authors investigated the evolution of 
the topological charge Q during the cooling process. They found that the lattice 
equilibrium configurations are dominated by configurations with Q = 0. However, 
if these configurations are cooled down, configurations with larger Q-values emerge. 
This demonstrates that the fluctuations on the lattice (at realistic values of the 
coupling strength) are large enough to change the topological charge Q. 

In this paper, we will elaborate a detailed picture of the response of an instanton 
to the presence of static charges on a "microscopic" level. We will thereby provide 
new information to clarify the role of the instantons for the dual Meissner effect 
rather than to estimate the importance of the (semi-) classical vacuum for the true 
ground state of Yang-Mills theory. 



3 The Meissner effect of solid state physics 



In 1933 Meissner and Ochsenfeld discovered that super-conductors expel magnetic 
flux. The microscopic ingredient which leads to this effect is the condensation of 
electric charge, which is carried by pairs of electrons [p9 |. For later comparison with 



the Yang-Mills case, we briefly review this effect by considering a Landau-Ginzburg 
theory pO] which is described in terms of the Lagrangian 



C = 1(0, + tA»(x)) <P(x)\ 2 + - Ux)A,(x) + V(<P 2 ) . , (1) 

where is the electro-magnetic field strength built from the gauge potential A^, 
and e is the electric charge. V(<p 2 ) is the tree-level Higgs-potential, which is minimal 
for a non- vanishing value </> 2 = §\ providing a condensation of the scalar field, jn(x) 
is the external current. The Lagrangian is invariant under U(l) gauge transforma- 
tions. The scalar field <fi(x) transforms homogeneously under these transformations 
and therefore carries electric charge. 

In order to mimic the scenario of a super-conductor, we assume that the scalar 
field forms a homogeneous condensate 0o of electric charge, which breaks the U(l) 
gauge symmetry. The classical equation of motion can be easily obtained from ([[]), 
i.e. 

d^u(x) - 2e 2 <gA v [x) = -j„(x). (2) 

In the gauge d^A^x) = 0, this equation can be transformed into an Helmholtz 
equation for the magnetic field £>j. In the absence of external currents, this equation 
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Figure 1: The Meissner effect: mag- 
netic field lines with a super-conducting 
cylinder present (solid lines) and with- 
out the cylinder (dashed lines). 



becomes 

d 2 B l (x) - 2e 2 (j) 2 B i (x) = . (3) 

This is the key equation to understand the Meissner effect. Eq.(|3|) tells us that the 
magnetic field exponentially decreases inside a super-conductor with slope given by 
the strength of the condensate 0o- Figure 1 qualitatively shows the behavior of the 
magnetic field lines, if a super-conducting cylinder is present. The main observation 
is that the magnetic field lines are expelled from the super-conducting region. 

Let us compare London's equation @ with the classical equation of motion of an 
SU(2) Yang-Mills theory. This equation can be derived from the Euclidean action 
which is defined in the appendix, i.e. 

d»F« u {A](x) - e abc F^[A)Al(x) = -g 2 f u (x) , (4) 

where e abc is the anti-symmetric tensor in 3 dimensions, and j„(x) are the external 
currents. The field strength, i.e. 

F;jA] = d,Al - d v Al + e ahc A\A c v , (5) 

is constructed from the non-Abelian gauge field A a {x). The external source j° is 
a real function of space-time and can be understood as the current generated by 
Euclidean quark fields (see (p3|)). We will adopt the Lorentz gauge, d^A®(x) = 0, 
and will study the impact of this current on instantons throughout this paper. 



6 



For a vanishing external current (J* = 0), eq.(|J) allows for non-trivial solutions, 
which are known as instantons [JTIJ. In the gauge d^A^x) = 0, their gauge field and 
the corresponding field strength is given by 

ainst a 2 ainst a V 

where 77" are the anti-symmetric 't Hooft symbols, i.e. 

•ft = i& = e°*, i,* = 1...3, (7) 

and p is the radius of the instanton. Instantons therefore correspond to spots of 
non-vanishing field strength. 

Let us now briefly discuss the linear response of the gauge field, i.e. a®(x), to 
the external source j%(x). Decomposing A a {x) = AV' nst (x) + a^(x), we assume that 
the current is sufficiently weak implying that one can expand the e.o.m. to linear 
order in a a Ax). If we confine us to the region close to the center of the instanton, 
i.e. y/x* p, the e.o.m. becomes 



d 2 a a v (x) - \e ahc F^ 



A mst (x = 0) a b Jx) = -g 2 f u (x);. 



Since the matrix e F° possesses negative eigenvalues, eq. (|) describes screening 



of the field strength as well as ant i- screening depending on the color index a under 
consideration. The field strength at the instanton center sets the scale of the 
(anti-) screening length and therefore plays a similar role as the scalar field 0q in (0). 
We learn from eq. (H) that Meissner type effects occur due to the non-linear nature 
of the Yang-Mills e.o.m., if instantons are present in the ground state of SU(2) Yang- 
Mills theory. The crucial difference is, however, that these effects lead to a repulsion 
or an attraction of the color-electric field lines depending on the orientation of the 
instanton. 

In this paper, we will further investigate these effects resorting to numerical 
methods in order to go beyond the linear response theory. We will provide evidence 
that the electric flux produced by the charges is expelled from the instanton core, 
if we average over the instanton orientation. In this case, the above leading order 
effect obtained by the linear response approach drops out. 
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4 The dual Meissner effect of SU(2) Yang-Mills 
theory 



We will numerically investigate the solutions of the classical equation of motion 
(|J) (we set g 2 = 1 for simplicity) of SU(2) Yang-Mills theory, which is a non- 
linear second order partial differential equation. To this aim, we discretize the 
4-dimensional Euclidean space-time on a grid consisting of 31 4 points, and replace 
derivatives by the corresponding differences, e.g. 



where e M is the unit vector in //-direction, and where h is the grid spacing. The 
discretized partial differential equation is then solved by iteration. For defmiteness, 
we have chosen Neumann boundary conditions. We have checked that the impact of 
the boundaries on the field configurations under consideration is small (see below). 

The approach should not be confused with the lattice version of Yang-Mills 
theory [IJ. The grid is only a mathematical tool to solve the partial differential 
equation. For our purposes, we need not worry about local gauge invariance, but 
confine ourself to the definite gauge choice (^|). 

In order to get a first idea of the error due to the discretization (and to check 
the set up of the program), we cast the analytically known instanton configuration 



4.1 The numerical approach 



In the following, we will work in the Lorentz gauge 
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Figure 3: The color-electric field E\{r = x{) of an instanton for 22,3,4 = 0: The 
numerical solution (triangles) in comparison with the analytic result (H). 



A^ mst (x) (|B|) onto the grid, and checked to what extent the discretized version of 
the l.h.s. of (fD reproduces the zero at the r.h.s. (note that the external source j° is 
set to zero at the moment). Choosing 22,3,4 = (this selects a line which passes the 
center of the instanton) and e.g. a = 2 and v = 3, the deviation from zero is shown 
in figure 2 as function of x\ for several values of the grid spacing h. The result should 
be compared with the intrinsic scale which is set by the maximum field strength of 
the instanton (4 in units of the instanton radius). The most important observation 
is that the result becomes significantly better, if h is decreased. This confirms that 
the discretization works correctly. 



4.2 The instanton on the grid 

In this subsection, we study the grid version of an instanton configuration, which 
emerges as a solution of the discretized e.o.m. (|j) at zero external source (j^(x) = 0). 
The SU(2) instanton possesses eight zero-modes, since the instanton configuration 
breaks symmetries of the Lagrangian. Four translational, three rotational and one 
dilatational degrees of freedom do not lead to a change of the action. In the pseudo- 
particle description of the instanton, these zero-modes give rise to collective co- 
ordinates, which must be fixed in order to arrive at a definite configuration. For 
this purpose, it is convenient to employ a procedure which minimizes the action by 
the method of steepest decent. Then zero-modes do not contribute to the gradient 
configuration which is added in each iteration step to the actual profile in order to 
further reduce the action. This implies that the starting configuration of the iter- 
ation process completely fixes the collective coordinates of the instanton. Here we 
have chosen the analytic profile (Bp as the starting configuration. 
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Figure 4: The color-electric field E l in the x|/-plane produced by two color-electric 
charges in a state with one instanton present. The dashed lines indicate the field 
lines of the charges in an empty vacuum. 
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After convergence was obtained, we have calculated the color-electric field of our 
grid configuration with the help of a discretized version of the field strength tensor 
@. Choosing the color direction a = 1, we learn from (||) and (0) that the only 
non-vanishing color-electric field points in e\ direction. We then have compared the 
numerical result for the space-time dependence of this color-electric field, i.e. E\(x), 
with the analytic form (^|). The result in shown in figure ^. For h = 0.5 p, the 
field strength in the center of the instanton is somewhat underestimated, whereas 
the result is satisfactory for h = 0.2 p. However, the influence of the boundaries 
becomes visible at the edges of the grid in the latter case. 



4.3 Electric charges in empty space 



In this subsection, we study the color-electric field produced by static color sources 
in non-Abelian SU(2) gauge theory. For this purpose, we have chosen the external 
current at the the r.h.s. of the equation of motion (HI) to be 



fJx) = 8 al 8^ [^(f-fO - 8®{x-x 2 ) 



'ID 



This choice corresponds to two static color-electric charges of color direction a = 1 
which are located at the spatial positions X\ and x 2 . 

In order to get a feeling of the color-electric flux produced by these sources, we 
calculated the color-electric field E a=l using an empty space as starting configuration 
of the iteration process, which solves the partial differential equation (|]). To be 
specific, we have chosen 

x, = (4p,-4p,0), x 2 = (-4p,-4p,0), (12) 

where p is an arbitrary length scale at the moment and will be the instanton radius 
later. If only sources for one color direction is present, the non- linear terms of the 
Yang-Mills equation of motion (^ drop out, and the e.o.m. essentially reduces to 
a Maxwell equation in the relevant color channel. For a later comparison with the 
case with an instanton present, we numerically calculated the color-electric field E l 
in the :q/-plane at each grid point. The color-electric field lines show the behavior 
anticipated from the analogy to classical electrodynamics. 



4.4 Electric charges in the instanton background 



The picture changes drastically, if the vacuum contains an instanton configuration. 
We here only consider sufficiently weak external charges which do not completely 
deform the instanton and which therefore preserve the pseudo-particle character of 
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the instanton. In this limit, we are interested in the interplay of the color-electric 
field strengths of the external charges and the instanton as function of the collective 
coordinates. Note that in the presence of the external charges, the collective coor- 
dinates of an instanton in empty space do not correspond to zero-modes anymore, 
but a particular choice of these parameters exists, if the action is minimal. One 
possibility to constrain the collective coordinates is to include the term 

xjd'x (F; v [Ae*\ - F«M}) 2 (13) 

to the action. The field A ex is thereby the gauge potential of a single instanton in 
empty space with a definite choice of origin, orientation and radius. The parameter 
A acts like a Lagrange multiplier. For sufficiently small values of A, the motion of 
the instanton along the (former) zero-mode directions due to the external source is 
blocked, and the collective coordinates are dictated by the background field. Us- 
ing several (small) values of A, one verifies that the field strengths configuration, 
produced by the instanton and the charges, remains (almost) unchanged by the 
additional term (p~3|) . 

In order to be specific, the iteration process which minimizes the action was 
started with the analytic instanton configuration (|6|) (which is also used as the 
background field in (|i~3D) casted onto the grid. The upper picture of figure |] shows 
the result for t — 0, where the instanton develops its maximum strength. The color- 
electric field lines of the instanton, which point from the right to the left, are clearly 
visible. The length of the vectors on the grid points also give a rough impression 
of the instanton field strength profile (|6]). Some color-electric field lines are also 
shown (solid lines). For comparison, the color-electric field lines of the case without 
the instanton are shown as dashed lines, too. A field line is specified by the angle 
between the field line and the line connecting the two charges. This makes it possible 
to compare the situation with and without instanton. Our main observation is that 
the color-electric field lines of the charges are clearly expelled from the region where 
the instanton is located. 

Let us study this effect, if the location of the two charges is changed to 

x x = (-4p,-4p,0), x 2 = (-4p,4p,0). (14) 

The result is shown in the lower picture of figure |j. The upper charge is the source 
of the color-electric field, the lower charge is the drain. One observes that the color- 
electric field lines which spread out the upper charge are pushed back from the 
instanton, whereas the lines of the drain charge are pulled in. The net effect of the 
instanton is to align the external electric field lines with its internal orientation. 

One might argue that the deformation of the color-electric field lines of the 
charges in the presence of the instanton is solely due to the superposition of the 
instanton field and the field from the charges. Figure 5 shows that this is not the 
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Figure 5: The color-electric field i? 1 
in the xy-plsme of the charge instanton 
system with the color-electric field of 
the unperturbed instanton subtracted. 
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case. In this picture, we have subtracted the color-electric field of a pure instanton 
(without any charges present) from the full color-electric field. If the total field is 
simply a superposition, one should recover the Maxwell type field lines (this is what 
would happen in classical electrodynamics). The result of this subtraction is com- 
pletely different from the field distribution expected from classical electrodynamics. 
This is possible due to the non-linear nature of the classical equation of motion 

5 Random instanton orientation 



In the previous section, we have investigated the distribution of the electric field 
strength, produced by a charge anti-charge pair, for a definite choice of the instanton 
orientation. No evidence for the dual Meissner effect was found so far, since the 
external color-electric field lines are attracted or pulled back from the instanton 
core depending on the orientation. Here, we will study whether a net repulsion 
of the color-electric field lines takes place, if the we average over the instanton 
orientation. For this purpose, we will calculate the space-time distribution the field 
strength squared F^ U F^ U , where we average over the orientation of the instanton and 
keep its position fixed. The case of the random instanton orientation is interesting 
from a physical point of view for the following reason: a realistic description of 
the Yang-Mills ground state resorts to a liquid of instantons [TE\. In this liquid, 
the orientation of a particular instanton is arbitrary. The small amounts of action 
which are necessary to move the instanton along its former (i.e. single instanton) 
zero-mode directions are provided by the entropy of the medium. 

Let us first study the dependence of the instanton interaction with the external 
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Figure 6: The instanton interaction Si 
as function of the instanton orientation 
9. 



charges on the orientation O of the instanton, i.e. 



X 2 + p 2 ' 



O 



cos 9 — sin 6 
sin 9 cos 9 




(15) 



where O ab is a particular orthogonal matrix parameterized by the angle 9. Since 
we will switch on electric charges in color a = 1 direction, it is sufficient for our 
purposes to study matrices O ab which rotate the instanton color-electric field in the 
£i£2-plane only. The instanton interaction Si is defined by 



Si — Sic — Si 



inst 



s, 



charges j 



(16) 



where Sic is the action of the instanton-charge system, and S inst and Surges is the 
action of a single instanton and the charges, respectively. If the charges are weak 
and far located from the instanton origin, the interaction Si can be calculated in an 
i.e. 



elegant way [16 



S"i oc r?^ v 



a ext. 



(17) 

where F°£ xt " is the field strength of the external charges at the instanton origin, 
and where the superscript indicates that it was assumed that the this field strength 
is slowly varying over the instanton region. If we align the charges in ei-direction 
symmetric to the X2-axis, their color-electric field at the instanton center is also 
oriented in e\ direction. This implies that S® oc cos# for the particular choice of 
orientation ([15]). The numerical result for Si is shown in figure 6. The strengths 
of the charges actually yield S 'charges I ' Si nst = 6.7%. The distance of one charge to 
the instanton center was 1.7 p. A fit according Si ~ a cos 9 + b is also shown in 
figure H A cos ^-modulation of the data, as predicted by (0), is clearly visible. In 
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Fi gure 7: (-F° l/ -F^) 1 ' 2 (xi, X2, Xo,3 — 0) for several values of Xi (cases a-c) of the 
charges in empty space (solid) and of the instanton-charge system with (F 2 ) 1 / 2 of 
an unperturbed instanton subtracted (dashed). In the case of the instanton-charge 
system, we have averaged over the instanton orientations. 



addition, one observes a non-vanishing offset b towards negative values. This offset 
can be understood as follows: if we put the field configurations of the charges and 
of the unperturbed instanton on the grid as starting configuration, the method of 
steepest decent further reduces the action in order to finally arrive at a solution of 
the Yang-Mills equations of motion. Since action of the starting configuration does 
not depend of the instanton orientation, the angle average of Si is negative. From 
this argument, it is obvious that the offset b cannot obtained in the linear response 
approach, which is the basis of eq. (|T7|). 

Finally, we will average over the #-angle, which defines the instanton orientation 
with respect to the charges. We will study the variation of the field strength squared, 
i.e. F" (x), in space-time rather than the integrated quantity Si. Since F 2 of a 
single instanton decreases like 1/r 8 for large distances r from the instanton center 
and since F 2 of a single charge asymptotically behaves like 1/r 4 , it is more convenient 
for illustration purposes to study the quantity \[F 2 . In the present investigation, 
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we have located the charges at 



%i = (P,-P,0) 



X 2 



-p, -p, o) 



'It 



Figure 7 shows y F 2 of the charge anti-charge pair (solid lines) as function of X\ 
(^3,0 — 0) for several values of x 2 , i.e. case a x 2 = — 0.4 p, case b x 2 = — 0.2 p and 
case ci2 = 0. This result is compared with yfF 2 ^ produced by the instanton charge 
configuration (dashed lines), where y/F 2 of a single instanton was subtracted and 
where over the orientations of the instanton was averaged. The latter quantity is 
therefore a measure of the field strength which exists on top of that of an unperturbed 
instanton. One clearly observes that y/F 2 produced by the external charges is 
suppressed at the space time region occupied by the instanton. In order to quantify 
this result, we introduce the ratio 



k{x 2 ) 



2 peak 
charges \ ^2) 



peak 



2 peak 
charges \ 2) 



(19) 



where \fW' Charges an d VF 2 ^ are the maximum values of y/F 2 of the charges 
in empty space and of the subtracted charge instanton system, respectively, at a 
given line (x 2 ,x 3)0 = 0). The ratio k directly measures the suppression of the 
external color-electric field inside the instanton. The numerical calculation reveals 
the following result 



x 2 





-0.2 


-0.4 


-0.6 


K 


92% 


84% 


63% 


28% 



One finds that a net suppression of the external field inside the instanton occurs, if 
over the instanton orientation is averaged. The suppression increases towards the 
instanton center. 



6 Discussions and conclusions 



Among many other things, Polikarpov and Veselov studied the instanton medium 
which emerges in lattice Yang-Mills theory, if an equilibrium configuration is cooled 



down [27]. They found in particular, that this cooled medium exhibits a non- 



vanishing string tension, which, however, is only 5% of the full string tension. 
Whether this residual confinement is due to some excitations still present on top 
of the instantons at the cooled lattice or due to an intrinsic property of the instan- 
ton medium is not clear yet. 
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In order to provide new information concerning the interesting question whether 
the instantons play a role for the dual Meissner effect, we have studied the impact 
of an instanton on the color-electric flux produced by two static, point-like color- 
electric charges. The classical equation of motion of a SU(2) gauge theory was 
numerically solved in Lorentz gauge. The color-electric field lines of the charges were 
investigated in detail in an "empty" vacuum as well as in a state where an instanton 
is present. We found that the external color electric field lines are expelled or 
attracted depending on the orientation of the instanton relative to the charges. The 
net effect is that the instanton aligns the external field lines which is in accordance 
with the result of the leading order of the linear response theory ||16|| . However, if 
we average over the orientations of the instanton, this leading order effect drops out. 
In sub-leading order, a substantial repulsion of the color-electric field lines from the 
instanton core region survives. 

The numerical outcome can be qualitatively understood as follows: in an "empty" 
space, the electric flux spreading out the charges distributes over the whole space in 
order to avoid large concentrations of field strength which amounts for a large action. 
If only charges of unique color are present, the scenario is the same for the Abelian 
(classical electro-dynamics) and the non-Abelian case. However, the non-linearity 
of the Yang-Mills equation of motion (in fact the topological structure of the the- 
ory) allows for non-trivial configurations, i.e. instantons, which are local minima of 
the Yang-Mills action. These configuration must exhibit a definite strength of the 
color-electric and color- magnetic field (see (^)) in order to minimize the action. Any 
distortion of the instanton field strength, e.g. produced by the presence of external 
charges, yields a sudden increase of the action. This is the reason that the color- 
electric field is expelled from the instanton core region, once the dominant dipole 
interaction of the instanton with the charges drops out due to a random instanton 
orientation. 

Whether this effect is a pre-cursor of quark confinement giving rise to a non- 
vanishing string tension, if the Yang-Mills ground state is modeled as a strongly 
correlated medium of randomly oriented instantons, is an interesting question, which 
is left to future studies. 
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Appendix: Euclidean gauge theory 

The Euclidean partition function of a SU(2) gauge theory, fermions included, is 
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described by the functional integral 



Z = JVq Vq ] VA^ exp j - J d l 



x 



—F; v F; u + q\i<]) + im-Alt a ) q 



(20) 



where the gauge fields A a ^ are considered as real fields, and where the fermion fields 
q transform under the fundamental representation of the SU(2) gauge group which 
is spanned by the generators t a . g is the gauge coupling constant and m is current 
quark mass. The action is invariant, if the fields transform under a global 0(4) 
symmetry, i.e. 

A a ; = A, V A: , q' = S(A)q , (21) 

which is the pendant to the Lorentz symmetry in Minkowski space. The orthogonal 
matrices A^ and the matrices 5* (A) are spanned by the anti-symmetric generators 



Ufu,, e.g. A^v = [exp{— The matrices S'(A) satisfy 

5(A) 7M 5 t (A) = A, ulu , 



(22) 



where 7 M are the hermitian Dirac matrices which fulfill {7^,7^} = 25^. From (p2[), 
it is obvious that 

j;(x) = q\x) 1 ,t a q(x) (23) 

is a hermitian current which transforms like a 0(4) vector. The zeroth component of 
this vector is the charge density which is the subject of investigations in this paper. 

Finally, let us check that the coupling of the current (|) to the gauge fields A a ^ 
is correctly chosen, i.e. a gauge transformation does not induce imaginary parts in 
A ^. One easily observes that the transformation 



q'(x) = U(x)q(x) 



U{x) e 5^(2) 



leaves the action invariant, if the gauge fields transform according 
t a A a '(x) = U(x)t a A a (x)U ] (x) - iU(x)d^(x) . 



(24) 



(25) 



It is easy to verify that the transformation (p5|) leaves the gauge field A ^ real. The 
contribution of the interaction of the quark current and the gauge fields to the 
Euclidean action is therefore 



Si, 



int 



- d*x j a Jx)A a Jx) . 



(26) 



It is straightforward to derive the equation of motion (^) from the action in (|20|) by 
taking the functional derivative with respect to the gauge fields A%(x). 

We here provide the right hand side of the equation of motion (Q), i.e. the 
current, and calculate the gauge fields. First information on the gauge fields can be 
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obtained by taking the divergence of the equation of motion. Using the e.o.m., a 
direct calculation yields 

d*fi{x) = e^A^x^x) . (27) 
In this paper, we study the response of the gauge fields to static charges, i.e. 

J» = fa'Eti (28) 

i 

Eq. fl2"T|) then tells us that the zeroth component of the gauge fields at the position 
of the charges either vanishes, i.e. A£(t, aT £ ) = 0, or that the gauge fields are oriented 
parallel to the charge color vector, i.e. A^(t, Xi) oc pf. 
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